Fast readout of a single Cooper-pair box using its quantum capacitance 
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We have fabricated a single Cooper-pair box (SCB) together with an on-chip lumped element 
resonator. By utilizing the quantum capacitance of the SCB, its state can be read out by detecting 
the phase of a radio-frequency (rf) signal reflected off the resonator. The resonator was optimized 
for fast readout. By studying quasiparticle tunneling events in the SCB, we have characterized the 
performance of the readout and found that we can perform a single shot parity measurement in 
approximately 50 ns. This is an order of magnitude faster than previously reported measurements. 
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I. INTRODUCTION 

Superconducting devices based on Josephson junc- 
tions have successfully been used in many different kinds 
of applications, including very sensitive magnetometers 
based on superconducting interference devices (SQUIDs), 
bolometric detectors, mixers, and parametric amplifiers. 
They have also been suggested to be strong candidates 
as building blocks for a quantum computer—"—. They are 
easily fabricated with standard lithographic techniques 
and can be integrated with other electrical circuits. This 
gives them the potential to be scalable. Devices utiliz- 
ing charging effects include very small tunnel junctions. 
Such Coulomb blockade devices, are also widely used in 
measurements, for example, the single electron transis- 
tors (SET)^. The radio frequency version of the SET is 
the worlds most sensitive electrometer^. 

The single Cooper-pair box (SCB)^"— is one of the 
simplest Coulomb blockade devices, involving a single 
Josephson junction. SCBs are very sensitive to the pres- 
ence of quasiparticles which suggested its use as potential 
radiation detector—. The presence of quasiparticles can 
be measured by detecting the charge on the SCB island 
using an external SETi^. In this paper, we characterize 
an intrinsic method for reading out the SCB which re- 
lies on the curvature of its energy bands. This methods 
is both faster and is predicted to have less backaction 
then using a SET. The curvature of the energy bands of 
the SCB (with respect to gate charge) gives rise to the 
so called quantum capacitanc o^^i^^ and has been utilized 
in a number of experiments, for example in the measure- 
ments on longitudinal dresses states of a driven SCBiS,ii£. 
It has also been used to study the ground state of two 
coupled qubit^ii, and to study quasiparticle poisoning 
of a SCB qubilji^. In Ref. ^^ they used a resonator 
with a bandwidth of 200 kHz which limited the speed 
the their measurements. Here we use the random tun- 
neling of quasiparticles to characterize the performance 
of the quantum capacitance readout. We show that we 
can measure the state of the SCB an order of magnitude 
faster than has previously been reported, including mea- 
surements using RF-SETs. We also show that we can 
prepare the SCB in a certain parity state with a high 



probability. 

This paper is structured as follows: in section [Hi we 
present the theory behind the readout technique of using 
the quantum capacitance. In section IIIIl we show how 
the samples were designed and fabricated and then in 
section ITVl we present the measurements done to charac- 
terize the readout. 



II. THEORY 

A. Cooper-pair box 

The sample under investigation is a SCB and is shown 
in Fig. [Tfa) & (b). The SCB consist of a superconduct- 
ing island connected to a large reservoir by a Josephson 
junction. The Josephson junction is made in a SQUID- 
configuration to allow the Josephson energy, Ej, to be 
tuned by applying a magnetic field through the loop. 
The SCB is also characterized by the electrostatic energy. 
Eel = EQ(n — UgY, which is the energy required to add 
an extra Cooper pair to the island. Here n is the number 
of Cooper pairs that have tunneled onto the island and 
Ug — CgVg/2e \8 thc normalized gate voltage, where Cg 
is the capacitance between the voltage source and the is- 
land. If the capacitance of the island, Cs, is small enough 
the Cooper-pair charging energy, Eq — (2e)^/2Ci;, will 
dominate over the Josephson energy, E'j, and the tem- 
perature, ksT. In this case the charge fluctuation on the 
island will be small. The number of excess Cooper pairs 
on the island, n, is then a good quantum number and the 
charge of the island can be well controlled by the external 
gate voltage, Vg. For Ej ^ Eq and < rt^ < 1, only two 
charge states will be of interest: |0) and correspond- 
ing to zero (n = 0) or one {n — \) extra Cooper pair on 
the island. Then the Hamiltonian of the Cooper-pair box 
can be written 

H = -]^EQ{l-2ng)cj,-^Ejn,, (1) 

where ax^ Uz are the Pauli spin matrices. Here we have 
ignored all state independent terms of the Hamiltonian. 
The two eigenenergies for this Hamiltonian are plotted 
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FIG. 1: (Color online) (a) A scanning electron micrograph of 
the resonator for sample A. The inductor can be seen spiral- 
ing around the plates of the capacitors with the upper part 
counter wound compared to the bottom, (b) A scanning elec- 
tron micrograph of the SCB for sample A. The SCB consists 
of a 5 nm long and 100 nm wide superconducting island. The 
SCB for sample B is the same except that the island is 8 
fim long in order to increase the gate capacitance. The is- 
land is connected to a reservoir through two small Josephson 
junctions in a SQUID geometry. The potential of the island 
is controlled by three different capacitive gates. One large 
rf gate, Cg^ , (on the top) is connected to the resonator and 
used for the readout. In addition, there is a dc gate, Cg'^, used 
to bias the SCB at the working point and a microwave gate, 
Cg™ , used for spectroscopy, (c) A schematic of the device. 
The parallel resonator consists of two metal layers separated 
by a thin insulating layer. The bottom layer is of Nb (within 
the dashed line), forming the inductor and bottom plates of 
the capacitors. Next is a 200 nm layer of silicon nitride (the 
dashed line) covering the whole sample and, finally, the top 
layer is Al (outside the dashed line) making the top plates of 
the capacitors as well as the SCB. (d) The energy of the two 
lowest energy eigenstates of the SCB, as well as the expecta- 
tion value of the (excess) charge on the SCB island for each 
state as a function of the normalized gate voltage, Ug. (e) 
The quantum capacitance, Cq, for the two eigenstates as a 
function of the normalized gate voltage, Ug. 



in Fig. [T]^d) as a function of Ug. At the degeneracy 
point Ug — 0.5, where the electrostatic energies of the 
two charge states cross, we get an avoided level crossing 
with a sphtting between the ground and excited state 



equal to Ej. In the same graph, we have also plotted the 
expectation value of the island charge (n) for each energy 
eigenstate. 



B. Quantum Capacitance 

At the degeneracy point the energy difference between 
the ground and excited state is, to first order, inde- 
pendent of the gate charge rig, which makes this the 
ideal bias point when the SCB is used as a qubit. At 
this point, the longest dephasing times are obtainecU^. 
Therefore, this is often called the optimal point. If we 
want to detect the state of the SCB sitting at the op- 
timal point, we cannot measure the charge, since the 
charges of the ground and exited state are the same at 
this point (see Fig. [IJd)). Although the charges are the 
same for the two states, the derivatives of the charges 
with respect to the gate voltage differ and can be used 
for readout. We can define an effective capacitance 
of the SCB by calculating Ceff = d{Qg) /dVg, where 
{Qg) = CgCjVg/Cs + 2e (n) Cg/C^ is the average value 
of the injected charge on the gate capacitor. This ef- 
fective capacitance, Ceff — Cgeom + C'q, will have two 
contributions, a geometric part Cgeom = CgCj / (Cg + Cj) 
consisting of the gate capacitance in series with the junc- 
tion capacitance, and a state dependent part that we call 
the quantum capacitance^^. This quantum capacitance, 
Cq, takes the following form 

^ C^{a^ + {l-2ng)^f^ C^ AE^ 

for the ground (g) and first excited (e) state, where a = 

Ej/Eq and AE = ^ Ej + E'^{1 - 2ngY is the energy 

difference between the two states at a given Ug. Cq is 
equal in magnitude but has opposite signs for the ground 
and excited state, with Cq being negative for the ground 
state. In Fig.[TJe), Cq is plotted against the gate charge, 
Ug, for the ground and exited state of the SCB. 

If we embed the SCB in a resonant tank circuit (see 
Fig.[IJc)) the effect of the quantum capacitance will be to 
shift the resonance frequency. The reflection coefficient, 
r = |r|exp(0), of the resonator has a constant magni- 
tude, i.e. |r| — 1, since there are no dissipative element 
in the resonator. The phase </>, however has a sharp fre- 
quency dependence and close to the resonance frequency 
of the tank circuit, = 27r/o — 1/ yjL{C + Cc), the 
phase, (/), can be approximated by the expression: 

= -TT - 2 arctan ( 2Q ^ ~ ] , (3) 
V ^0 J 

where uj = 27r/ is the angular probe frequency and Q — 
(C -I- Cc)/C^ZoOJo is the external quality factor of the 
resonator. Now, if the capacitance of the SCB is changed 
by a magnitude of ACq, the phase of the rf signal at 
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frequency /o reflected off the resonator will change by: 

A</.«-2arctan(^Q^^^^ . (4) 

Here we have treated the resonator purely classically. A 
full quantum treatment, including the SCB, resonator 
and transmission line can be found in the work of Jo- 
hansson et al.'^^ . It is shown that this method of readout 
is quantum limited^ meaning that no information is lost 
during the readout (no extra dephasing). 



C. Quasiparticles 

Quasiparticles are single-particle excitations of the su- 
perconducting condensate. Quasiparticle fluctuations 
have been studied as a source of noise in superconduct- 
ing devices for some time. Thermodynamic fluctuations 
in the quasiparticle number, also know as generation- 
recombination noise, is an important source of noise 
at intermediate temperatures^i. Time-resolved mea- 
surements of these fluctuations have shown very good 
agreement between theory and experimen t^^i^^ . At very 
low temperatures, were thermal quasiparticles should be 
suppressed, a significant population of quasiparticles is 
still observed in most experiments. The origin of these 
nonequilibrium quasiparticles is still unknown. However, 
it is clear that they remain an important source of noise. 

The most significant source of quasiparticle noise in 
a SCB is commonly referred to as quasiparticle poison- 
ing. When a nonequilibrium quasiparticle in the reservoir 
tunnels onto the SCB island, it shifts the potential of the 
island by ±e/Cs. The random tunneling of quasiparti- 
cles on and off of the island, therefore leads to a large- 
amplitude telegraph noise in the island potential. Quasi- 
particle poisoning has been extensively studie d^^i^^'^^ , in- 
cluding the observation of individual tunneling events in 
real tim o^^i^^ . It has been shown that poisoning is well 
described by a simple kinetic model, starting from the 
assumption of a finite density of nonequilibrium quasi- 
particles in the leads. In this model, the tunneling rates 
are then dictated by the relative energies of the even (no 
quasiparticles on the island) and odd (one quasiparti- 
cle on the island) states. If the energy of the odd state 
is lower than the even state, the quasiparticle can be 
trapped on the island (see Fig. |3l[a)). The energy dif- 
ference of the even and odd state is maximum at the 
charge degeneracy point, Ug = 0.5, where it takes the 
value 5E ^ Eq/A~ Ej/2 - A.; + A^. Here A^ and A^ 
are the superconducting energy gap of the island and the 
reservoir, respectively. The average occupation of the 
two states is then simply determined by the Boltzmann 
factor exp{-SE/kT). 



III. DEVICE DESIGN AND FABRICATION 



When designing the readout circuit there is a trade off 
between having a good signal (large phase shift) and a 
low backaction. The phase shift is roughly proportional 
to QCq/{C + Cc) (see equation |4]), meaning that for a 
fixed Q- value {i.e. measurement bandwidth) we want 
a small total capacitance in the resonator (and a large 
inductance). However, the voltage noise of the environ- 
ment will induce charge fiuctuations on the SCB, with a 
magnitude that scales with the same pref actor. Thus, by 
decreasing the total capacitance in the resonator, C-|-Cc, 
you can increase the sensitivity (larger phase shift) at the 
cost of a larger backaction on the SCB. 

The Q-value of the resonator sets the bandwidth, i.e. 
an upper limit on how fast we can measure. We designed 
the resonator to have an external Q-value of 100 which 
corresponds to a time constant of 50 ns at /o = 650 MHz. 
The internal Q-value is usually substantially larger and 
can therefore be ignored. 

The devices were fabricated in a multilayer process. 
Starting from a high-resistivity silicon wafer with a na- 
tive oxide, the wafer was first cleaned using rf back sput- 
tering directly after which a 60 nm thick layer of niobium 
was sputtered. To pattern the niobium, we used a 20 nm 
thick Al mask made by e-beam lithography and e-beam 
evaporation. The niobium was then etched in a CF4 
plasma (with a small flow of oxygen) to form the inductor 
and bottom plates of the capacitors (see Fig. [TJc)). The 
choice of niobium for the bottom layer made it possible to 
test the resonator in liquid helium (even with a top layer 
made of normal metal, e.g. gold). The Al mask was re- 
moved with a wet-etch solution based on phosphoric acid. 
Before depositing the insulator we cleaned the wafer in 
a 2% HF solution for 30 sec in order to remove most of 
the niobium oxide which has been found to degrade the 
Q-value of the niobium resonators. Using PE-CVD, we 
then deposited an insulating layer of 200 nm of silicon 
nitride. The silicon nitride layer covers the whole wafer; 
connections to the niobium layer are only made through 
capacitors. We chose silicon nitride since it is known to 
have low dielectric losses^. After using a combination 
of DUV photolithography to define bonding pads along 
with e-beam lithography to define quasiparticle trapsSS., a 
3/80/10 nm thick layer of Ti/Au/Pd was deposited by e- 
beam evaporation. Finally, the layer containing the SCB 
was made by e-beam lithography and two-angle shadow 
evaporation of lO-f-30 nm of aluminum, with 6 min of 
oxidation at 4 mbar. The thickness of the island (10 nm) 
was chosen to be much thinner than for the reservoir 
(30 nm) in order to enhance the superconducting gap, 
Ai > Ar, of the island compared to the reservoir. This 
was done to reduce quasiparticle poisoning. 
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IV. MEASUREMENTS 

A. Measurement setup 

The device was cooled in a dilution refrigerator with a 
base temperature of about 20 mK. For the readout, we 
used an Aeroflex 3020 signal generator to produce the 
rf signal. The signal was heavily attenuated and filtered 
and was fed to the tank circuit via a Pamtech circulator 
positioned at the mixing chamber. The reflected signal 
was amplified by a Quinstar amplifier at 4K with a nom- 
inal noise temperature of f K. The in-phase and quadra- 
ture components of the signal were finally measured us- 
ing an Aeroflex 3030 vector digitizer. A schematic of the 
measurement setup can be found in Fig. [51 



TABLE I: The extracted parameters for the two samples. 





Sample A 


Sample B 


/o 


676 MHz 


663 MHz 


Q 


128 


130 


L 


151 nH 


324 nH 


C 


251 £F 


97 fF 


Cc 


116 fF 


81 fF 


CI' 


0.2 fF 


0.3 fF 


Eq 


62 GHz 


48 GHz 


Ej 


7.2 GHz 


7.4 GHz 



B. Device characterization 




300 K 



4.2 K 



1.7 K 



0.6 K 



20 mK 



FIG. 2: (Color online) The measurement setup used in the 
experiment. The sample (inside the dashed line) is mounted 
inside a copper box on the mixing chamber of a dilution refrig- 
erator, which reaches a temperature of 20 mK. The readout 
pulse is created using a Aeroflex 3020 signal generator (with 
built in IQ modulation). It is heavily attenuated and cou- 
pled through a circulator located at the mixing chamber and 
filtered through two RLC waveguide filters before reaching 
the resonator. On the way up, the signal is first amplified at 
4.2K with a Quinstar amplifier with noise temperature of 1 K 
and then amplified again at room temperature with a Mini- 
Circuit amplifier. Before the signal is finally digitized using 
an Aeroflex 3030 vector analyzer (with a 33 MHz bandwidth), 
it is first filtered using an RLC bandpass filter to reject image 



In the following sections we will show measurements of 
two different devices, referred to as sample A and sample 
B. In order to characterize the devices and extract param- 
eters, the resonator was first measured using a network 
analyzer. From the measured phase response we could 
extract the resonance frequency, /o, and the Q-value of 
the resonator. We could extract both the Cooper-pair 
charging energy, Eq, and the maximum Josephson en- 
ergy, Ey^"'^ by conventional spectroscopy, while applying 
a perpendicular magnetic field through the SQUID loop 
of the SCB, and thereby tuning Ej. Finally by measur- 
ing Cq as a function of gate charge, we can extract the 
rf gate capacitance and the total resonator capacitance. 
From the expression for the Q-value and the resonance 
frequency, we can then extract the parameters for the 
individual components of the resonator. The extracted 
parameters for the two devices are presented in Table HI 
The values of the parameters for sample B corresponds 
reasonable well to the geometrically identical device that 
was used in Ref. [2^. The gate capacitance for sample 
A roughly agrees with what you would expect from the 
5 /.tm island compared the to the capacitance of the 8 
/im island in sample B. If we the insert the parameters 
from Table U into Eq. (g]) we calculate an expected phase 
shift of 11 deg for sample A and 30 deg for sample B 
between the two parity states when biased at the degen- 
eracy point. 



C. Readout performance 

We have fabricated and measured a number of samples. 
All devices we have measured so far have been poisoned, 
meaning that quasiparticles cause switching between the 
two parity states of having an even or odd number of 
quasiparticles on the island. The switching between these 
parities happens on the time scale of a few microseconds. 
Although this is far from ideal for many applications, it 
has given us a way to characterize the readout. When 
a quasiparticle tunnels on or off the island the quantum 
capacitance of the SCB will change, and can thus be de- 
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have also performed pulsed measurements of the state. 
We send down a Gaussian pulse, with a length defined as 
the full width at half maximum (FWHM), and measure 
the phase of the reflected pulse. In order to optimize 
the signal to noise of the measured response, we used 
a so-called matched filter, where the time traces of the 
measured in-phase and quadrature component (I and Q) 
are multiplied with a Gaussian template (with the same 
shape as the one generated by the signal generator). The 
product is then integrated and a single value for is ex- 
tracted. This is done for both I and Q and we can cal- 
culate the phase hy (f) = arctan(Q//). We perform 100 
000 of these measurements and make histograms of the 
measured phase and see two peaks centered at different 
phases corresponding to the two parity states. We fit 
a double Gaussian (both with the same standard devia- 
tion) to the histograms (see Fig. |4lja,b)). We define the 
signal-to-noise ratio (SNR) of the measurement as the 
peak separation divided by the standard deviation. This 
is performed for different pulse lengths and we extract the 
SNR as a function of the measurement time, Fig. |3Uc,d). 
The SNR roughly follows the expected square root de- 
pendence on the measurement time. To reach a SNR 
> 1 we need a pulse length of the order of 50-100 ns 
in both samples. The shortest measurement time was 
here limited by the time constant of resonator which was 
about 60 ns for both of the samples. 



FIG. 3: (Color online) (a) The energy as a function of Ug for 
the even (solid) and odd (dashed) parity of the SCB island. 
The odd state energy is shifted by the difference in the su- 
perconducting gaps of the island and the leads, Ai — A;. If 
Eq /4 — Ej/2 > Ai — A; the island will form a trap for quasi- 
particles at rig = 0.5. (b) The phase shift as a function of the 
gate charge, rig, for the even (solid) and odd (dashed) state. 
We also show the average measured phase for sample B, which 
shows contributions from both parity states, while repetitively 
sweeping the dc gate at a repetition rate of 5 kHz. The rate 
of the dc sweep was fast compared to the quasiparticle tun- 
neling rate, such that the quasiparticle tunneling probability 
during one period was low. The SCB then spent most of the 
time in the even state, which has on average a lower energy 
than the odd state, (c) A typical time trace of the measured 
phase, from an identical device to sample B, when sitting at 
the even degeneracy point, Ug = 0.5. The SCB spends most 
of the time in the odd state {(f) — deg) but makes short 
excursions to the even state {(j) = —30 deg). 



tected as a change in the reflected phase from the res- 
onator. The measured time averaged phase for sample 
B is shown in Fig. El^b) as function of Ug. In Fig. [3jc) 
we show the time dependence of the phase measured at 
the degeneracy point, Ug — 0.5. Most of the time the 
SCB is in the odd state with a phase shift of about 
degrees, however now and then the extra quasiparticle 
escapes the island and the SCB spends short periods in 
the even state with a phase shift close to -30 degrees. 
While sitting at the degeneracy point, Ug = 0.5, we 



D. Quasiparticle relaxation and state preparation 

We know that by sitting at the even degeneracy point 
the system will eventually relax into the odd parity state 
given that Eq/A -Ej/2 > A, - A/ (see Fig. O. We 
wanted to study how fast this process is and to what ex- 
tent you can prepare the SCB in a certain parity state. 
We start by letting the system equilibrate biased at the 
odd degeneracy point, i.e. Ug = 0, thereby preparing the 
even state. We then pulse the gate to the even degeneracy 
point (at t = 100 /is) and observe the dynamics. From 
a long time trace, including 1000 pulses, we divide each 
repetition into 0.5 /is increments. We extract the average 
phase from each increment and each repetition. We then 
make a histogram of the phase for each increment as a 
function of time (see Fig.[5l^b)). We fit the histograms to 
a double Gaussian and extract the occupation probabil- 
ity of the even and odd state as a function of time (see 
Fig. EIc)). Sitting a the odd degeneracy point, Ug — 0, 
there is an equilibrium probability of more than 80% of 
being in the even state and when we pulse to the even 
state most of the probability is preserved. Eventually the 
system equilibrated and then the probability is reduced 
to 20%. From the relaxation of the probabilities we ex- 
tract an equilibration time of 2.8 jis. As a comparison, 
we show the average time trace (see Fig.[5Ua)) where we 
have taken the average of the full time traces from dif- 
ferent pulses. From this we extract a relaxation rate of 
3 /zs in good agreement. Since the quasiparticle relax- 
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FIG. 4: (Color online) Sitting at Ug = 0.5 we perform pulsed 
measurements and extract the phase of each pulse. For three 
different pulse lengths we repeat the measurements many 
times and make histograms of the phases (bars). In (a), we 
show histograms for measurements on sample A and in (b) for 
sample B. Sample A was designed to have a larger total ca- 
pacitance for the resonator and a smaller rf-gate capacitance 
than in sample B. This reduces the backaction of the measure- 
ment at the expense of lowering the SNR. To the histograms, 
we fit a double Gaussian (solid line) with the same standard 
deviation for the two peaks. We then define the SNR of the 
measurement as the peak separation divided by the standard 
deviation. The above procedure is repeated for different pulse 
lengths, ranging from 20 ns to 1 /xs. The SNR is then plotted 
as function of the pulse length for sample A in (c) and for 
sample B in (d). We see that the SNR roughly follows the 
expected square root dependence on the measurement time. 
Since sample A was designed to have a lower backaction, we 
expect it to have a lower signal. The calculated phase shift 
is 11 deg for sample A and 30 deg for sample B. This value 
is however calculated for very low probe powers. In these 
measurements, the phase shift is reduced by the relative large 
measurement power. 
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FIG. 5: (Color online) Measurement of the quasiparticle re- 
laxation in sample B. (a) The average of 1000 traces of the 
measured phase after pulsing the gate from the odd degener- 
acy point {rig = 0) (at time 100 fis) to the even degeneracy 
point {ug = 0.5) and continuously monitoring the phase as a 
function of time, (b) Histograms of the measured state as a 
function of time. The measured phase for each repetition is 
divided into 0.5 ^s increments. We then make a histogram 
of the average phase for each increment, measured from the 
start of the pulse. This produces a histogram as a function 
of time, (c) Occupation probability of the even (lower) and 
odd (upper) state as function of time extracted from the his- 
tograms in b. After pulsing to the even degeneracy point, the 
probabilities are inverted suggesting that we should be able 
to prepare the system in the even state with high probability 
even if the device has significant quasiparticle poisoning. 



CONCLUSIONS 



ation rate is much longer than the operation and readout 
times this suggests that, even if the device is poisoned, 
we should with high probability be able to prepare the 
box in the even state and perform useful measurements. 



We have fabricated and tested two samples with a SCB 
together with an on-chip lumped element resonator. The 
resonators were optimized for readout speed, with a Q- 
value around 100. We have characterized the readout by 
employing the effect of quasiparticle poisoning and found 
that for readout pulses of length 50-100 ns we get a SNR 
greater then 1. 
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